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Abstract 

Doubled topological phases introduced by Kitaev, Levin and Wen sup- 
ported on two dimensional lattices are Hamiltonian versions of three dimen- 
sional topological quantum field theories described by the Turaev-Viro state 
sum models. We introduce the latter with an emphasis on obtaining them 
from theories in the continuum. Equivalence of the previous models in the 
ground state are shown in case of the honeycomb lattice and the gauge group 
being a finite group by means of the well-known duality transformation be- 
tween the group algebra and the spin network basis of lattice gauge theory. 
An analysis of the ribbon operators describing excitations in both types of 
models and the three dimensional geometrical interpretation are given. 

1 Introduction 

Topological quantum fields theories (TQFT) in three dimensions describe a vari- 
ety of physical and toy models in many areas of modern physics. The absence of 
local degrees of freedom is a great simplification, it often leads to complete solv- 
ability [TJ |2]. Perhaps the most recent territory, where they appeared to describe 
real physical systems, is that of topological phases of matter, being e.g. responsible 
for the fractional quantum Hall effect [3j. Since the idea of fault-tolerant quantum 
computation appeared in the literature pE], TQFT's are also important in quan- 
tum information theory. These new applications also enhanced the mathematical 
research, led to classification of the simplest models p]. 

Due to their topological nature, TQFT's admit discretization yet remaining an 
exact description of the theory given by an action functional on a continuous mani- 
fold. One large class thereof are the so-called BF theories, whose Lagrangian density 
is given by the wedge product of a (d — 2) form B and the curvature 2-form F of 
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a gauge field [6]. We will deal here with a special class of three dimensional theo- 
ries, which describe doubled topological phases and restrict our attention to discrete 
gauge groups G. The context they appeared in, in recent physics literature [H [7] , are 
Hilbert spaces of states in two dimensions and dynamics therein, which are bound- 
ary Hilbert spaces H of the relevant TQFT's. Operators acting on H correspond 
to three dimensional amplitudes on the thickened surface. In this paper we will ex- 
plain this correspondence, which was proved for the ground state projection recently 
[H], and provide the geometric interpretation of the ribbon operators, which create 
quasi-particle excitations from the ground state. This is a step towards extending 
the correspondence to identify the ribbon operators as invariants of manifolds with 
coloured links embedded in them in the TQFT. 

The emergence of topological phases from a description of microscopic degrees 
of freedom are modelled by the lattice models of Kitaev j3] and Levin and Wen [7j. 
Since they generically have degenerate ground states and quasi-particle excitations 
insensitive to local disturbances, they are also investigated in the theory of quantum 
computation [9], their continuum limit being closely related to the spin network 
simulator [TUJ [TT] . The ground states were extensively studied in the literature, 
their MERA (multi-scale entanglement renormalisation ansatz) [121 13J and tensor 
network representations [Hj have been constructed to study e.g. their entanglement 
properties pj)J ITS) . Finding the explicit root of these structures in lattice gauge 
theory and TQFT can help to understand their physical properties. 

Lattice gauge theories admit seemingly very different descriptions. A state can be 
represented by assigning elements of the gauge group to edges of the lattice. The dual 
description in terms of spin network states, where edges are labelled by irreducible 
representations (irreps) of the gauge group and vertices by invariant intertwiners are 
also well known since [17] . To name an application, this description turned out to 
provide a convenient basis for most approaches to modern quantum gravity theories 
[T5], US]. In this paper we will show in detail how these dual descriptions give rise to 
Kitaev's quantum double models in one hand and the spin net models of Levin and 
Wen on the other. Then the ribbon operators in both models and their identification 
will be discussed. 

The organisation of the paper is as follows. In the next section, we introduce 
the Turaev-Viro models via the example of BF theories. In section 3, we briefly 
introduce the string net models of Levin and Wen on the honeycomb lattice V in the 
surface S and recall the proof [5j that the ground state projection is given by the 
Turaev-Viro amplitude on S x [0, 1]. The boundary triangulations of S are given by 
the dual graphs of V decorated by the labels inherited from the "initial" and "final" 
spin nets. In section 4 the duality between the states of the Kitaev model and 
the string nets will be shown by changing the basis from the group algebra to the 
Fourier one. By using this duality and an additional projection, we will obtain the 
electric constraint operators of the string net models. The matrix elements of the 
magnetic constraints are also recovered provided that the local rules of Levin and 
Wen hold. We explain that they do in all BF theories, which is a strong motivation 
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in their favour for the case when the gauge group is finite. In section 5, we discuss 
ribbon operators and give their three dimensional geometric interpretation in terms 
of framed links in the Turaev-Viro picture. Finally, a summary is given with a list 
of questions for future research. 

2 Turaev— Viro models 

In three dimensions both the F(A) (the field strength) and B fields of BF theory can 
be be considered to be forms valued in the Lie algebra of the gauge group G. The 
action can then be written as j M tr(BAF) with tr being an invariant non-degenerate 
bilinear form on the Lie algebra and M is a smooth, oriented, closed three-manifold. 
We may start from the case when G is a semi-simple Lie group relevant in particle 
physics theories and gravity, A being the connection in the principal G-bundle over 
M. In three dimensions the "space-time" separated form of the Lagrangian has 
the structure BjdAk + A DjAk + B Fjk(A) where j,k are spatial indiceaj. The 
first term is the standard kinetic term, the second implies the (Gauss) constraint of 
gauge invariance, the third the vanishing of the (two-dimensional) field strength (A 
and B are Lagrange multipliers), D and d stand for the covariant and the exterior 
derivatives, respectively. 

Since locally the solution of the constraints is given by a pure gauge A; = g~ x dig 
(g : S — > G smooth function, with S being the spatial hypersurface) one may 
discretize the theory by introducing a lattice on the spatial surface and quantize 
the remaining degrees of freedom: the holonomies (elements of G) describing the 
coordinate change between faces of the lattice. They correspond to the edges of the 
dual lattice, which is constructed by placing a vertex inside each face and connecting 
new vertices, which were put inside neighbouring faces. This dual lattice is the 
starting point of the models in [I] , the electric constraints are the remainders of the 
Gauss constraint and the magnetic ones are the remainders of the flatness constraint. 
For a detailed exposition see e.g. |20j. 

The partition function of the above BF theory is formally obtained by taking 
the functional integral over the fields A an B of the phase associated to the classical 
action 



It is not so easy to give this definition a precise sense, but for the moment it is 
not necessary to go into further details. What does matter is that there exists a 
consistent way to discretize the partition function by considering an oriented tri- 
angulation Ma of the manifold M by assigning two Lie algebra element Bi, Q t to 
each edge i in Ma. The generator B i can be thought as the integral of B along 
the edge i, whereas f2« as the logarithm of the group element corresponding to the 

1 There is not necessarily physical time in the theory, one can do this decomposition for Euclidean 
signature as well. 
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holonomy around the edge J^. Then the Feynman integral in (JTJ) can be replaced 
by Y\ i J dgi J dBie 1 $ trBlQi . The B integrals will yield Dirac deltas 8 gi> x and one can 
now proceed with decompositions in terms of irreps of the gauge group. This way 
one ends up with a discrete state sum instead of the original Feynman integral, 
where each state is the triangulation coloured with irreps and its weight is given by 
the precise final form of the amplitude (examples are given below). The structure 
of the partition function (amplitude) for a prototypical theory, the Ponzano-Regge 
model [21] corresponding to G = SU(2), reads 

^ M A)=^n^ri( 6 ^' ( 2 ) 

where (6j) is the Wigner 6j symbol of SU(2) depending on the 6 irreps decorating 
the edges of the tetrahedron t, d^ is the dimension of the irreps ji assigned to the 
edge i and the sum ranges over all states, that is, all possible colourings of the edges 
with irreps. It turns out that this state sum is well defined and independent of the 
chosen triangulation for a large class of mo delfi For a systematic derivation of this 
state sum from action functionals, see [22] or section 2.3 of [23] . 

The Ponzano-Regge partition function (j2J) is formally independent of A, but is 
divergent. However, the Turaev-Viro (TV) model [24] . a regularized version thereof, 
has a well defined partition function, given by 

^lRll{i; * % 

] % t 

where the underlying algebraic structure is the quasi-triangular Hopf algebra SU q {2) 
with q = exp(27u//c), fceZ fixed, j 6 [0, 1, . . . , k — 1] denotes the irreps of SU q (2), 
dj G C is the so-called quantum dimension of j, the constant d is defined by 
d = J2 k dl, the quantity in the brackets is the quantum 6j symbol and V is the 
number of vertices of the triangulation. One finds the precise definitions of all the 
quantities along with the algebraic properties assuring consistency and triangulation 
independence [21] of the amplitudes in [21]. We will briefly mention the origin of 
the latter property in the next section. The final fact for this introductory section is 
about the form of the amplitude ([3]) for manifolds with non-empty boundary. The 
associated boundary triangulation, whose edges are decorated by labels {j'}, derives 
from a given triangulation in the 3d bulk and is kept fixed. The amplitude reads 

*tv[m\ {/};,] = y, "- v IK IK 1 IK n ft I £} • w 

j int i V v' t ' 

2 To be more precise, one needs to introduce the dual complex by putting vertices inside every 
tetrahedron, connecting those vertices which were put in neighbouring tetrahedra and a vertex 
should be singled out on the boundary of each dual face. Then the procedure to get is the 
following: Take the dual face corresponding to i. Multiply the holonomies along the boundary 
edges of this dual face starting from the vertex singled out in a circular direction determined by 
the orientation of i (say, by the right hand rule). The logarithm of this group clement is f^. 

3 This is one way to define a TQFT rigorously. 
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where V is the number of internal vertices, the index i ranges over internal edges, 
%' over boundary edges, v' over boundary vertices (each boundary vertex is the end- 
point of an internal edge, j v i is its colour), t over all tetrahedra and the summation 
is done for internal edge labels only, while those on the boundary are kept fixed. 

Note that there is a quasi-triangular Hopf algebra associated to finite groups as 
well, the so-called Drinfeld (quantum) double T>{G) [25J. There, the dimensions dj 
as well as the 6j symbols can be obtained from the representation theory of the 
group G. 



3 String nets 



Levin and Wen [7] started off from the algebraic structure underlying the above 
models (consistent set of 6j symbols and quantum dimensions), which serves as the 
algebraic data in defining TQFT's. Taking these data for granted, they constructed 
a two dimensional lattice model, which we will now introduce briefly. Consider a 
surface S with a fixed oriented honeycomb lattice F embedded in it. The Hilbert 
space is spanned by all possible decorations of the edges with labels j; we will refer 
to them as irreps (of SU q {2) or the finite group G) as we will not need to treat 
the most general TQFT's. The Hamiltonian is a sum of two families of mutually 
commuting constraint operators, 



H = - 



(5) 



where the first sum is over all vertices, the second is over all plaquettes of the lattice. 
A v = Nij k for k being the irreps decorating the edges adjacent to the vertex v 



(the numbers 



E N are referred to as fusion coefficients between the irreps: 
The magnetic constraints are written as a sum 



B n 



(6) 



over irreps and the action of the individual terms is 




ET?bg*h r?ch*i rpdi*j r?ej*k jpfk*l rpal*g 
s*h'g'* s*i'h'* s*j'i'* s*k'j'* s*l'k'* s*g'V* 

g'h'i'fk'V 




(7) 



while its action on the rest of the state supported on the honeycomb lattice F is 
trivial. The numbers F^ n are the 6j symbols, part of the algebraic data of a TQFT, 
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i* denotes the irreps dual to i. Changing the orientation of an edge is equivalent to 
changing its label i to its dual i*. Levin and Wen use a different normalisation than 
that of ©, ©: 

Before proceeding, let us write down an important algebraic property of the F 
symbols: 

N 

Erpmlq rpjip rpjs*n _ rpjip T?riq* /q\ 
kp* n runs* Ikr* q*kr* mis* \ ) 

n=0 

This identity is called the Biedernharn-Elliot identity or pentagon equation, which 
holds in every TQFT. In the concrete examples mentioned above, they can be proved 
by the definition of the 6j symbols as connecting the two different fusion channels 
of recoupling irreps (graphically encoded by ( J26l) in section 4.2)0. 



3.1 Reconstructing 3d geometry 

In our work [8] we recovered a three dimensional Turaev-Viro invariant [2EJ EZ] from 
the algebra of Levin and Wen. We associated geometric tetrahedra to the algebraic 
6j symbols, where the edges are decorated with irreps from the 6j symbols. In 
that a convention needs to be adopted, e.g. the upper row should correspond to a 
(triangular) face of the tetrahedron and labels in the same column should correspond 
to opposite edges. In the examples we are looking at there is always a normalization 
of the 6j's such that they have the same symmetry as the tetrahedron. Orientation 
of edges can also be taken care of in a consistent manner, we shall however omit 
them for most of what follows. Now we can translate the Biedernharn-Elliot identity 
to geometry. The following picture arises: 




where the two configurations (three tetrahedra joined at the edge n and two 
tetrahedra glued along the triangle (irq)) correspond to the left-hand and right- 
hand side of ([9]), respectively. This is a cornerstone in proving triangulation inde- 
pendence of the amplitudes ©dlD and shows that whenever tetrahedra are glued, 
labels corresponding to internal edges have to be summed over. 

Now, one can try to find the geometric counterpart of the operator (171). Con- 
structing the dual (triangle) graph T of the honeycomb lattice T such that a dual 
edge inherits the label of the original edge it corresponds to (recall that edges of 

4 by means of using two different ways of coupling five irreps. 
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the original and the dual graph are in 1-1 correspondence in two dimensions), we 
proved the equality [8]: 

<rf } | n^l r ?"' } > = z tv[s x [o, i],f?W>] • (io) 
p 

The lhs. means the matrix element of the operator FT B p between two spin nets, 
that is, the honeycomb lattice To = T decorated by labels {j} and another copy of 
r\ = T decorated by {j"}- The rhs. is the Turaev Viro amplitude (j3J) of the three 
dimensional manifold S x [0, 1] with fixed triangulations on the two boundaries 
given by the dual graph Tj with the labels inherited from r$. In figure a) below, the 
dashed lines show a part of f ^ . Let us concentrate on the middle vertex in the figure 
corresponding to the dual triangle abc. There is an F symbol corresponding to that 
vertex from all three operators B p of the three plaquettes sharing that vertex. To 
each F we associate a tetrahedron and they induce the internal triangulation of S x 
[0, 1] depending on the order how the B p operators are multiplied one after the other. 
These different orders of multiplication correspond to different decompositions of the 
prism (built from translating the triangle abc in To to the corresponding one a"b"c" 
in Ti) into three tetrahedra. The fact that they commute is nicely reflected by 
the independence of the TV amplitude on the internal triangulatiorjj. The figure 
b) shows a part of the amplitude corresponding to F^l^/F^^^F^,^,, . The labels 
with one bar a', V , d are to be summed over when the full amplitude is written in 
accordance with the fact that the underlying edges are the internal edges of the 
triangulation of S x [0, 1]. Summation of the labels comes from the sum in (J6j). 



a) b) 




Note that FT B p = FJ B p Y[ v A v if we naturally define the 6j symbols to be zero 
whenever a triple (ijk) corresponding to a dual triangle of a tetrahedron has Nijk = 
0. This way, what we found is that the TV amplitude gives the ground state projec- 
tion. For the precise matching of all the weights to those of (j3J) and the consistency 
of the full amplitude, see the last section of [8]. 

4 Duality and the quantum double lattice models 

Now, we shall restrict our attention to the case, when the TQFT is given by the 
structure of the double T>(G) of a finite group G and show the equivalence between 

5 Note that for proving this equality it is necessary that the coefficients of are given by ((6j). 
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the lattice models of Kitaev p[] based on that structure and the corresponding string 
net models. Our method relies on the duality in the underlying lattice gauge theory 
[T7] . Essentially the same idea was employed also in the very recent paper [28J. 

The Hilbert space of the Kitaev model (and that of a lattice gauge theory) 
is spanned by the group algebra basis {\gi, 92, • • • , <7e), 9i G G}, supported on an 
oriented lattice T with E edges V vertices and F faces. The scalar product is given 

(91,92, ■ ■ ■ , 9E\hi, h 2 ,..., h E ) = S gi!hl Sg 2ih2 . . . 5 gE: h E (11) 

The Hamiltonian consists of two families of sums of constraint operators, which are 
projections and mutually commuting. We will follow the strategy of imposing the 
electric constraints first and find the corresponding operators in the string net model 
of Levin and Wen. Then we will study the action of the magnetic constraints in the 
range of the set of electric constraints, and determine their matrix elements in the 
dual basis, recovering the magnetic operators in the string net model this way. 

The basic idea is the well-known expansion of any function / G L 2 (G) with G 
being any compact Lie group 

/(j)^^c:A( S ) (12) 

j m,n 

in terms of irreps j (D^ are the representation matrices and c J mn are coefficients). 
The statement is known as the Peter- Weyl theorem. We now define a new basis 

{\j1J2, ■ ■ ■ jE,oti,a 2 , ■ ■■0£ E ,(3 1 ,P 2 , ■ ■ -,Pe)} (13) 
by means of the scalar product 

(9i,92, ■ ■ ■,9e\3i,32, ■ ■ ■ ,jE,ati,a 2 , . ..a E ,p 1 ,p 2 , ■ ■ ■ , Pe) 



= D^( gi ) ai(}l D^(g 2 ) a ^ 2 . . . D^{g E ) aE p E 

The Qjj (Pi) denote the target (source) index of the oriented edge i and they range 
over the dimensions of the irreps j. We will need a linear combination of this basis 
defined in the following way. Consider all elements with fixed irreps ji,j 2 , ■ ■ ■ ,3e- 
For every vertex v of T take a three index-tensor I v , where the indices range over 
the dimension of the irreps associated to the three edges (the honeycomb lattice is 
trivalent) incident to v. Then contract all indices with the corresponding ones in 
(fT5|) . A simple example corresponding to the theta-graph is given in the figure. 



D Jl (gi) Pl . 
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For these states associated to the graph T we use the notation 



\ji,h, ■ ■ -,jE,Ii,h, • (14) 

4.1 The electric constraints 

Let us recall the electric constraints of the Kitaev model. They are written in terms 
of the following operators 

L 9 (i v) • I h v | ^ f | ■ • ■ 5"^ • • • ) if i points towards v, 

' » • • • / \ | . . . 1 . . . ) otherwise. 

The local gauge transformation acting at vertex v reads 

A g (v) = J[W(i,v) , (15) 

(the product is over edges incident to the vertex v) and the electric constraint is the 
projection defined as the average of the latter over the group: 

A (v) = t^Y.Mv)- (16) 
1 1 g&G 

Note that the range of the set of electric constraints are gauge invariant states, that 
is, they are invariant under Y[ v A 9v (v) with arbitrary tuple (gi, g 2 , . . . , gv) G G v , as 
shown by the following calculation: 

' ' feeG ' ' feec ' ' /i'eG 

Hence, at each vertex, the projection (Tl6|) implements gauge invariance. Let us see 
how this is done in the general set of states (THI) also called spin networks. The 
action of a gauge transformation A v (g) on a spin network can be determined by 
rewriting the scalar product as 

(9i, 92,-.., g E \ E A9 (v)s,s>\S') = (S\A9(vy\ 9l ,g 2 , ...,g E ). (17) 

S' 

Let us write down this action explicitly for a vertex v whose incident edges are 
oriented outwards and use a simpler notation \g ly g 2 , g$) for a generic state supported 
on T with 1, 2, 3 being the labels of the edges incident to v. Let us also use a similar 
abbreviation j 2 , j$, I v ) for \S) as the remaining parts are not important for the 
case at hand. Since A 9 (vY = Yliev L 9 (i, vY = Yliev L g (i*, v) with i* denoting the 
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opposite orientation for the edge i we can write 

01, 02,03 

= (ji,h,33, Iv\ggi,g92,99a) = E I v ia2a3DJl (ggi)^i3 1 D 32 {gg2)a 2 2 D : > 3 (gg 3 ) a3 p 3 ■ ■ ■ 

01 ,02,03 

= E ( E ir 2a3 D j Kg)a 1J1 D j %g) a2l2 D^(g) a3l3 Y 

71,72,73 \ai ,02,03 / 

D^( gi ) JlPl D^(g 2 ) r2P2 D^(g 3 )^ 3 • • • = £ I^ 273 D^( gi ) J10l D^(g 2 )^ 2 D^(g 3 ) y3(33 ■ ■ ■ 

71,72,73 

In the above the dots • • • stand for the remaining part of the spin network, which 
is not affected. In the fourth equality the group homomorphism property of the 
matrices (D(gh) = D(g)D{h)) were used. The last equality is the definition of I£ 
as the quantity in the big parenthesis. 

The transformation rule of I v — > 1% means that I v G ji (g> j 2 (g) js- The Clebsch- 
Gordan series i ® j = Y2k Nijk k shows that (after a suitable unitary transformation) 
I v is block diagonal with j x ® j 2 ® J3 = Ylu ^jij2k^kj 3 i I and each block transforms 
according to an irreps / of G. The case I® = I v , g G G for all g G G corresponds 
to the trivial representation, which appears in the block decomposition iff Nj 1 j 2 j :j 7^ 
0. We see now that gauge invariance at vertex v is achieved by acting with the 
projection that projects into the invariant subspace of the decomposition. This 
should corresponds to the projection of Levin and Wen. Assuming iV^ < 2 for 
every triple of irreps ensures that there is one unique gauge invariant tensor / (for 
coupling three irreps), an intertwiner, so that the invariant subspace of Y[ v A 9 V V is 
spanned by 

{bij J2> • • • iJe) '■ Nj it j kt j t = 1 for all (2, k, I) incident to a vertex} 

and I is understood to be at every vertex contracting all a, (3 indices of ( fT3l) . For 
these state associated to the honeycomb graph Y and only for these, we will use the 
notation IS'). 

A shorter way to arrive at invariant spin network states is to consider a generic 
gauge invariant state supported on Y in the group algebra basis. These are the 
so-called cylindrical functions \l/ G L 2 {G E ) with the invariance property 

<*|0i, 92, ■ ■ ■ , 9e) = V(gi, 92, ■ ■ ■ 9e) = 

(18) 

^(hit^gthisx)- 1 , h(t 2 )g 2 h(s 2 )-\ h(t E )g E h(s E )- 1 ) 

for every (hi, h 2 , . . . , h v ) G G v where t(i) (s(i)) denotes the target (source) vertex of 
the edge i. It can be shown that the spin network states constitute an orthonormal 
basis in this Hilbert space |29j. 
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4.2 The magnetic constraints 

To recall the construction of the magnetic operators of the Kitaev model, we define 
auxiliary operators associated to pairs (i,p) where p is a face (plaquette) of T and 
% G dp is an edge on the boundary of p: 

T 9 (j,p) : |0i,02, ■■■i9e) >-> 5 g ±i g .\gi, g 2 , ...,gs) , (19) 

where we have g (g^ 1 ) in the argument of the Dirac delta, when p is to the right 
(left) of the edge i oriented forward. The magnetic constraint is the special case 
g = 1 of the operator 

6 

B M = E W Thm ^P). (20) 

To adapt to the string net model we took T to be the honeycomb lattice. After 
straightforward calculation one finds the action of -Bi(p) to be given by 

\9i,92, ■ ■ ■ ,9e) i-> $9 pm g m -9p(6),i\9i, 92, ■ ■ ■ ,9e) , (21) 

whenever all edges p(l),p(2), . . . ,p(6) bounding the hexagonal face consecutively 
point to the counterclockwise direction. Should a boundary edge p{l) point to the 
opposite direction, g p ^ needs to be replaced by g~,* in the above expression. To 
proceed we write down the Plancherel decomposition of the Dirac delta function, 
which reads 

<W.. S e,i = $>itr(I>W • . .fl,)) = J2 d MD j (gi)D j (g 2 ) . . . D*(g 6 )) ■ (22) 

3 j 

Each term in the above sum is the scalar product of a spin network based on 
a two-valent graph, the hexagon with \gi, g 2 , . . . , g§). This is a spin network of 
only bivalent vertices^! "evaluated" on the same group elements that appear in 
the bounding plaquette p of the spin network \S). So we may write the action 
of B p : \S) i — > dj\S,p, j) , where the state \S,p,j) is a generalized spin network 
with double lines inside the plaquette p. The notion, used also in [7], nonetheless, 



still requires proper definition. Were the use of the local rules 

*( — )= *( -Or ) (23) 

$ ( = ( ^ (24) 

)= M( T^>T ) (25) 

*( ^ )= £„*r*( 33 ) (26) 



6 The only intertwiner for bivalent vertex is the trivial Dirac delta connecting identical repre- 
sentations. 
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of [7j allowed, we could just refer to the calculation given by formula (CI) in that 
article, which gives the expansion of \S,p,j) in terms of bona fide spin networks 
(S\gi, #2j • • • j 9e)- We could then just take it as the definition and we would be 
done. However, to argue in favour of these local rules in the Kitaev model, we 
need to get back to the theory in the continuum. It has been mentioned that in 
the case when the group G is a Lie group, the electric constraints are the lattice 
versions of the Gauss constraint that imposes local gauge invariance. This was 
explicitly justified in the previous section. Turning to the flatness constraint: any 
flat connection has trivial holonomy g 7 [A] = Vexp(J^ A) along a closed curve 7 that 
is contractible (otherwise we could contract the curve to the point, whose curvature 
would be proportional to the generator of the holonomy). The converse is also true, 
to every decoration of T with group elements satisfying the constraint (l2"Tj) for all 
plaquettes, there exist smooth flat connection(s) in the manifold T is embedded into. 
Suppose that we have constructed one for the embedding surface of the honeycomb 
lattice. Then a spin network state 



with any graph T' makes sense and it is invariant of the homotopy class of the graph 
T'. This justifies ( 1231) . The connection is flat, so the holonomy along a contractible 
curve is 1, TrD^l) = dj, which gives ( 1241) . There is no non-trivial intertwiner 
between two different irreps, whereas the lhs. of ( 1251) is a composition of invariant 
maps with i — > j included, so that rule also holds. Finally, we can smoothly contract 
the edge with label m in (1261) without changing the value of (1271) and then we have 



where the middle equality is a property of intertwiners and the rightmost formula 
coincides with the inner part of the rhs. of f[2"B"]) . when its middle edge with label n 
is contracted. Note that we have omitted also the representation matrices for the 
irreps i, j, k, I as they are not affected by the above, as well as the other parts of the 
spin networks. 

Let us summarize what we have achieved. If we have a Lie group G and impose 
gauge invariance on the honeycomb lattice T, the matrix elements of the magnetic 
operators in the Kitaev model in the spin network basis {\S)} can be done in two 
step. First, one constructs a smooth connection in the manifold in which T is 
embedded. Then uses the local rules for transforming the spin network (1271) as in 
formula (CI) of [7J. During this process, the group elements also change as we 
deform the edges, whose holonomies are these group elements, but in the end, we 
can deform all edges to their original location. This way we find a linear combination 




(27) 




(28) 
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of the spin network states {\S)} corresponding to the magnetic constraint given by 
the expression ((7j). 

Nevertheless for finite groups the local rules are, even if well motivated, postu- 
lates. The magnetic operator has been derived in a more direct way by introducing 
some auxiliary degrees of freedom in the very recent paper [28] . 



5 Ribbon operators 

In section 3.1 we have been studying the ground state, the constraints that it sta- 
bilizes and the projection from the Hilbert space into the ground state as a three 
dimensional TV amplitude. One of the main physical interests, however, is the 
string-like excitations, the ribbon operators, which correspond to quasi-particles. 
We are going to sketch the corresponding preliminary results to illustrate that the 
logic which, worked for the ground state projection, provides us with the three 
dimensional interpretation of these quantities as well. 

A general ribbon in the spin net model is a string running along a certain path 
in the honeycomb lattice. The corresponding operator has the following structure 



where k runs through the vertices of the string, is the label of the third edge 
adjacent to the fc-th vertex, which is not part of the string. The label s« is the 
"type" of the string. The index structure of the 6j symbols are given by 

F e *l k k ~,l , if P turns left at Ik; 
F s k = { "fr 1 * (30) 

F ., i 1 , if P turns right at I k ; 

s Vfc-i 

where Ik is the /c'th vertex of the string. The Q matrices in a string operator have 
the index structure 




' v ik v. 



v., ~ if p turns ri s ht ' left at I k> : k+i; 



^ = <j ^ 4fc , if P turns left, right at I k , I k+1 ; (31) 



S SkSk+1 ■ Id, otherwise , 

and generically Q l j kl are matrices (so they have two more indices, which are sup- 
pressed above). We would like to proceed as in section 3 and find a TV amplitude 
that a string operator describes. Before asking what the Q matrices correspond to, 
let us see what geometry we find by passing the description to the dual graph and 
gluing a tetrahedon whenever there is an F symbol. 
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In the left figure we have depicted a part of a string, indicating the dual graph 
along. In the following we will mean this line when referring to the string, and 
we will mean the collection of dual triangles (shown by green dashed lines in the 
left figure) when referring to the ribbon. In the middle figure we took the ribbon 
and drew a tetrahedron over each triangle it consists of, as dictated by ( 1291) . The 
decoration of the edges follow the index structure of the operator. The edges with 
the same label belong to the same edge of the spin net, so they are to be glued. This 
results in the right figure. 

We may interpret the above in the following way. There is the path ABCD . . . 
in the dual graph fg , which is a continuous line of dual edges of the ribbon that 
correspond to edges of T, which connect vertices with different turning directions of 
the string (left and right). There is an analogous path A'B'C'D' . . . in f\ j ' } . The 
gluing dictated by the algebraic structure of the operator is such that the line in 
f f ^ winds around the one in f ^ exactly once during each segment of the line. 
For each such segment an Q matrix is present in the form of the operator and the 
notation Tj (i = 0, 1) refers to the initial and the final string net. 

The observables in the TV model are typically ribbon graphs, fat graphs or 
links embedded in a manifold, over the labels of which, there is no summation in 
the amplitude [26, 27J. They are invariant under isotopy transformations. This 
property is ensured by the precise form of the braiding matrices, which then satisfy 
the Yang-Baxter equation. The latter equation seems to be related to eqn. (22) 
of [7] in the spin nets. However, the precise relation and the identification of the 
braiding matrix in the TV models with an expression of Q as e.g. the work [30] 
suggests should be found for a complete equivalence. 

5.1 Kitaev's ribbons 

The ribbon operators are present in the Kitaev model as well pi]. 
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A prototypical example shown in the figure is given by a strip between a path 
along the edges of the original lattice (thick lines) and a neighbouring path in the 
dual lattice (dashed lines). It can be composed of elementary operators associated 
to triangles, which connect sites, that is, pairs of a plaquette and a vertex on its 
boundary. In the figure sites are indicated by green dotted lines. One elementary 
building block (i,p) is a triangle, which is composed of an edge i and a dual vertex 
(which corresponds to a plaquette p). The other elementary building block (j,v) 
is also a triangle composed of a dual edge (which corresponds to the original edge 
j) and a vertex v. The associated operators depend on elements of the double 
T>(G), which can be represented by pairs of group elements (g, h). The two types of 
elementary ribbon operators read 

W^ 9 \i,v) = 5 9 , x L%v) W'^ 9 \j,p) = T°-\j,l) . (32) 

Recall that the lattice is assumed to be oriented so these formulae make sense. 
The composition of these elementary operators into a long ribbon is done by the 
comultiplication, which is given by 

ur(h,g) _ ST^ w (hi,ai) w (h 2 , g2 ) (Kg) w ithu; (M) -5 5u h S i 

hi,gi 

(33) 

It is desirable to express these ribbons in the spin network basis to recover their 
corresponding matrix elements in the spin net model. However, there are several 
obstacles, which should be overcome to accomplish this task. In [7j, there are ad- 
ditional local rules rules to reduce a generalized spin net containing ribbons, to the 
basis {|S')}, see the beginning of section 4. In order for this to work, one should find 
a generalized spin network representation of the above operators. Another difficulty 
comes about when the dual string crosses the original one. In this case, the ele- 
mentary triangles overlap and the corresponding comultiplication operations do not 
commute. One needs to find a consistent rule to define their comultiplication in a 
non-ambiguous way. Note that the simplest ribbon operator in the spin net model, 
that is the one which winds around one hexagon, is easily found to correspond to 
(1201) . We find the following equality: 

B a (p) = Y,tr(D J (9))Bi > . (34) 

3 

The procedure to get it is doing the comultiplication for the six elementary opera- 
tors all corresponding to the second type in ( |32l) to arrive at 5 gi g2 g 3g4 g5g6 i- Here 
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the Qi are the group elements corresponding to the edges in the group algebra basis 
{\9ii92i ■ ■ ■ i9e)}- Then one draws a generalized spin network representation corre- 
sponding to the Plancherel decomposition of the Dirac delta as shown in the figure 
(similarly to those for the magnetic constraints) and resolves it to the spin network 
basis by using the local moves. It is however not straightforward to generalize itQ. 




6 Summary and outlook 

In this paper we have been studying the lattice models of Levin, Wen and Kitaev 
from two perspectives. On the one hand we identified the ground states and the 
constraint operators of these models in case the underlying lattice is the honeycomb 
and the gauge group is a finite group. This has been achieved by changing the basis 
from that of the group algebra, that is, when edges are decorated by group elements, 
to the Fourier basis. This basis is spanned by the matrix elements of the irreps. A 
special linear combination by means of invariant intertwiners at the vertices has been 
shown to provide the range of all electric constraints and the projection at individual 
vertices has been identified with the projection to the invariant subspace. Then, the 
magnetic operators in the group algebra basis have been shown to correspond to 
those in the spin net model once the local rules postulated in the latter are satisfied. 
We gave an argument in favour of them from lattice gauge theory with continuous 
gauge group. 

A second focus of the paper was on mapping the spin net to the Turaev Viro state 
sum. We have used the idea of building up simplicial manifolds by tetrahedra with 
edges decorated with irreps corresponding to 6j symbols in the algebraic expressions 
of operators in the spin net model. This provided the three dimensional geometric 
interpretation for the ribbon operator. Having a precise TV amplitude identified 
with the ribbon operator in the spin net needs further investigation. 

One would also like to match these ribbon operators also in the model of Kitaev 
and the spin net of Levin and Wen. However, finding generalized spin network 
representations of the previous so that one could reduce them to the spin network 
basis is not straightforward. 

In a series of papers (3TJ [32l [33] , families of g-deformed 'spin network automata' 
were implemented for processing efficiently classes of computationally-hard prob- 
lems in geometric topology -in particular, approximate calculations of topological 

furthermore, the argument given in section 4.2 in favour of the local rules is also lost, since 
the underlying connection here is not flat. 
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invariants of links (collections of knots) and of closed 3-manifolds. A prominent 
role was played there by 'universal' unitary braiding operators associated with suit- 
able representations of the braid group in the tensor algebra of (SU(2) q ). Traces of 
matrices of these representations provide polynomial invariants of SU (2) q-colored 
links (actually framed links), while weighted sums of the latter give topological in- 
variants of 3-manifolds presented as complements of framed knots in the 3-sphere. 
These invariants are in turn recognized as partition functions and vacuum expecta- 
tion values of physical observables (Wilson loop operators) in 3-dimensional Chern- 
Simons- Witten (CSW) Topological Quantum Field Theory [TJ. As is well known 
(see e.g. [6], the review [31] and the original references therein), any 3d TQFT of 
BF-type can be presented as a "double" CSW model, on the one hand, and the 
square modulus of the Witten invariant for a closed oriented 3-manifold equals the 
TV invariant for the same manifold, on the other. 

The remarks above make it manifest that the efficient (approximate) quantum 
algorithms proposed in [31] [32], [33] could be extended in a quite straightforward way 
to the string- net ground states and ribbon-like excitations framed in the "naturally 
discretized" double SU (2) CSW environment given by the TV approach, as we have 
done in the present paper. Work is in progress in this direction. 
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